Draft version March 4, 2013 

Preprint typeset using MIjiX style emulateapj v. 5/2/1 1 



(N 

o 

(N 
bJQ' 

< 
00 



CO 
^' 

6 

> 

in 
cn 

o 



X 



PROGENITOR-EXPLOSION CONNECTION AND REMNANT BIRTH MASSES FOR NEUTRINO-DRIVEN 

SUPERNOVAE OF IRON-CORE PROGENITORS 

Marcella Ugliano 1 , Hans-Thomas Janka 1 , Andreas Marek 1 , and Almudena Arcones 2,3 

Draft version March 4, 2013 

ABSTRACT 

We perform hydrodynamic supernova simulations in spherical symmetry for over 100 single stars of solar 
metallicity to explore the progenitor-explosion and progenitor-remnant connections established by the neutrino- 
driven mechanism. We use an approximative treatment of neutrino transport and replace the high-density 
interior of the neutron star (NS) by an inner boundary condition based on an analytic proto-NS core-cooling 
model, whose free parameters are chosen such that explosion energy, nickel production, and energy release 
by the compact remnant of progenitors around 20 M are compatible with Supernova 1987 A. Thus we are 
able to simulate the accretion phase, initiation of the explosion, subsequent neutrino-driven wind phase for 
15-20 s, and the further evolution of the blast wave for hours to days until fallback is completed. Our results 
challenge long-standing paradigms. We find that remnant mass, launch time, and properties of the explosion 
depend strongly on the stellar structure and exhibit large variability even in narrow intervals of the progenitors' 
zero-age-main-sequence mass. While all progenitors with masses below ~15M© yield NSs, black hole (BH) 
as well as NS formation is possible for more massive stars, where partial loss of the hydrogen envelope leads to 
weak reverse shocks and weak fallback. Our NS baryonic masses of ~1.2-2.OM and BH masses >6 M Q are 
compatible with a possible lack of low-mass BHs in the empirical distribution. Neutrino heating accounts for 
SN energies between some 10 50 erg and ~2 x 10 51 erg, but can hardly explain more energetic explosions and 
nickel masses higher than 0. 1-0.2 M . These seem to require an alternative SN mechanism. 

Subject headings: supernovae: general — stars: neutron — stars: massive — stars: evolution 



1. INTRODUCTION 

One of the main goals of the long-standing quest for a bet- 
ter understanding of the supernova (SN) explosion mecha- 
nism is the establishment of a theoretical connection between 
the properties of progenitor stars and those of SNe and their 
remnants. Understanding the progenitor-remnant systemat- 
ics would not only mean a better definition of SNe as the 
end point of massive-star evolution, it could also serve as 
input for continuative theoretical studies like binary popula- 
tion synthesis calculations. Moreover, it would provide indis- 
pensable information and the theoretical basis for deducing 
constraints from observations on the physical processes that 
trigger and power the onset of the explosion at the center of 
the SN. In particular, the explosion mechanism plays a cru- 
cial role for determining the energetic and dynamical prop- 
erties of the SN blast wave. It thus controls the conditions 
for explosive nucleosynthesis and is of central importance for 
answering long-standing astrophysical questions: Which pro- 
genitors give birth to neutron stars (NSs), which ones to black 
holes (BHs)? How do SN properties depend on the progen- 
itor metallicity? Are the observed masses of compact rem- 
nants compatible with modeling results? Do observed SNe 
provide evidence for more than one explosion mechanism? In 
particular, considering neutrino heating as the most widely 
favor ed process to trigger the explosions (B ethe & Wilsonl 
119851) . the specific question arises, which range of SN en- 
ergies and nickel masses can be explained by the neutrino- 
driven mechanism. 
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A popular way to model the blast wave initiation for study- 
ing SN nucleosynthesis and determinin g compact remnant 
masses are piston-driven explosions (e.g. JWooslev & Weaver! 
1 1995b IZhang et"aT1 12008, and references therein), in which 
a piston is placed at a chosen position in the pre- or post- 
collapse stellar core and moved with a prescribed, time- 
dependent velocity to mimic the infall of the selected mass 
shell and the explosive expansion of the newly formed SN 
shock. Besides the initial mass cut — fallback of nonescap- 
ing material may shift the final mass cut — also the energy 
deposited by the piston is prescribed as a free parameter. 
Alternatively, internal energy "bombs" have been used, in 
which a defined amount of thermal energy is released by 
an artificial increase of t he temperature in a certain volume 
dAufderheide et alj|1991l) . 

Investigating spherically symmetric stellar core collapse 
by hy drodynamic simu l ations with simplified neutrino treat- 
ment, |Q' Connor &"Otl (1201 lb identified a single parameter, 
the "progenitor bounce compactness" £2.5 of the innermost 
2.5 M of a star, whose value above a certain threshold they 
considered to be indicative for an unlikely explosion and thus 
BH formation. Progenitors with £2.5 ^ 0.45 were found 
to reach the BH formation limit typically within <0.8s af- 
ter bounce because of their tremendous mass accretion rates. 
Such stars also turned out to require a particularly high time- 
averaged neutrino-heating efficiency to develop an explosion 
against the huge ram pressure of the infalling core material. 
Nevertheless, the exact value of £2.5 where the bifurcation of 
the behavior occurs, appears somewhat arbitrarily chosen and 
not supported by strong theoretical arguments. In particular, 
the BH formation limit could already take place at a lower 
v alue of ^2.5 ■ 

iBelczvnski et alj d2011l) and iFrver et al.1 (12012b made the- 
oretical predictions of the NS/BH mass distribution by ap- 
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plying a simple analytic concept for estimating the explo- 
sion time, explosion energy, and fallback mass on the basi s 
of structural properties o f the progenitor stars (iFrverl 120061) . 
iFryer & Kalogeral (1200 ll) derived theoretical BH mass distri- 
butions by employing the explosion energy as function of the 
progenitor mass from a fit to results of a small set of four 
2D models of neutrino-driven explosions (iFrverll 19991) . They 
considered different fallback masses in dependence on the 
(parametized) fraction of the energy released by the SN mech- 
anism that is used to unbind the outer layers of the progenitor 
stars. 

In this work we adopt a different approach. We perform hy- 
drodynamical simulations in spherical symmetry for a large 
set of 101 progenitor models, following core collapse and 
bounce, postbounce accretion, the possible onset of an ex- 
plosion and subsequent neutrino-driven wind phase of the 
cooling proto-neutron star (PNS), and the development of the 
SN blast wave over hours to days until the fallback mass 
is determined. In order to power the explosions, we refer 
to the neutrino-driven mechanism, whose viability seems to 
be supported by recent 2D si mulations with detaile d energy- 
dependent neut rino transport dMarek & Jankall2009l) and gen- 
eral relativity dMiiller et al.l l2012aUbh . although an ultimate 
robustness of the model predictions may require 3D simula- 
tions and may hinge on remaining uncertainties of the high- 
density neutrino opacities and equation of state (EoS). Since 
in our spherical simulations explosions by neutrino heating 
must be triggered artificially, we excise the high-density core 
of the nascent NS and replace it by a time-dependent, ana- 
lytic cooling model. This core model yields the time evolu- 
tion of the neutrino emission from the high-density interior 
of the PNS. Its free parameters are adjusted such that obser- 
vational characteristics of SN 1987A, i.e., the explosion en- 
ergy, E% 1/K ~ 1.3 x 10 51 erg, and ejected 56 Ni mass, Wnl87a ~ 
0.07 Mr, are reproduced for progenitor stars aro und 20 M Q 
(e.g.. iTanaka etaDfl999t lUtrobin & Chugaill201 ll) . This im- 
plies the assumption that SN 1987A was an ordinary case of a 
neutrino-driven explosion in which, e.g., exceptionally rapid 
core rotation did not play a crucial role for the collapse dy- 
namics. The NS core model with the same parameter values 
is then also applied to all other collapsing stars. 

Provided SN 1987 A was such a normal case and neutrino- 
energy deposition was responsible for its explosion, our ap- 
proach thus allows us to explore the response of other progen- 
itor stars to the input of neutrino energy released from their 
forming compact remnants, and in the case of successful ex- 
plosions to determine the blast wave and remnant properties 
of the stars. In the following we will give more information 
on our modeling strategy and tools, present our results, and 
discuss the implications. 

2. NUMERICAL SETUP AND INPUT 

We employ the explicit, finite-volum e, Eulerian multi-fluid 
code Prometheus dFryxell et al.ll 19891) to solve the ID hydro- 
dynamics equations on a spherical grid. We use typically 
1000 geometrically spaced radial zones with the outer bound- 
ary at a radius R Q \,. The latter is placed at a large distance (ini- 
tially at 150,000 km) that is not directly affected by the simu- 
lated stellar collapse and SN explosion. In particular it is suf- 
ficiently large to prevent the SN shock from leaving the com- 
putational domain during the simulated time. The advection 
of nuclear species is treate d by the Consistent Mult i-fluid Ad- 
vection (CMA) scheme of iPlewa & Miillej ( 119991) . and self- 
gravity of the stellar plasma includes corrections of the New- 



tonian gravitational p otential for general r elativ istic effects as 
describ ed in detail in IScheck e t al. ( 2006|) and lArcones et al.l 
(120071) . 

The inner core of the PNS with a mass of M c = 1.1 M Q 
and densities well above those of the neutrinospheric layer is 
excised and replaced by a point mass at the coordinate ori- 
gin. The shrinking of the PNS is mimicked by a retraction 
of the closed inner boundary at Rn, to gether with the radia l 
grid, applying the function defined in lArcones et al.l (120071) . 
For the simulations discussed in this paper we use the values 
#j b = #(1.1 M Q at f~10ms p.b.)~ 60km, R{ b = 20km, and 
fib = 0.4 s for the initial and final boundary radii and the expo- 
nential contraction timescale, respectively. The chosen values 
yield a boundary co ntraction similar to the "standard case" of 
IScheck et alj ( 120061) . Hydrostatic equilibrium is assumed at 
the boundary radius R,\, between excised core and accretion 
envelope of the nascent NS. 

The cooling of the PNS interior and corresponding bound- 
ary luminosities are des cribe d by the neutrino core-cooling 
model presented in Sect. 12.11 The temperatures of the spectra 
of inflowing neutrinos (assumed to have Fermi-Dirac shape) 
are set equal to the gas temperature 7\ in the innermost grid 
zone, T V (R\\,) = T- u consistent with the fact that neutrinos 
in highly optically thick regions are close to local thermal 
equilibrium. The spectral degeneracy parameters are chosen 
to be r] Ve — r\y e = 77v t = for electron neutrinos v e , elec- 
tron antineutrinos v e , and heavy-lepton neutrinos v x (= v M , 
V/j, v T , v T ), respectively. (Using the values for local chem- 
ical equilibrium for Tj Ve and rj„ e does not cause any signif- 
icant differences.) Neutrino transport and neutrino-matter 
intera ctions in the compu tational domain are approximated 
as in IScheck et al.l (2006) by radial integration of the one- 
dimensional (spherical), grey transport equations for neutrino 
number and energy, assuming the neutrino spectra to have 
Fermi-Dirac shape with a calculated local temperature T v . 

The high-density plasma is described b y the tabulated 
microp hysical equation of state (EoS) of iJanka & Miillerl 
(119961) including arbitarily degenerate and arbitarily relativis- 
tic electrons and positrons, photons, and four predefined 
nuclear species (neutrons, protons, alpha particles, and a 
representative Fe-group nucleus) in nuclear statistical equi- 
librium (NSE). At temperatures and densities below NSE 
the thermodynamics are treated by the Helmholtz-EoS of 
iTimmes & Swestyl (120001 ). 

The NS cooling and its neutrino emission are followed for 
8-25 s until the shock has crossed 10 10 cm. At that time the 
power of the neutrino-driven wind, which pushes the dense 
ejecta shell enclosed by the expanding SN shock and the 
wind-termination shock, has decayed to a dynamically in- 
significant level. The models are then mapped to a grid ex- 
tending to larger radii, and the inner boundary is moved from 
below the neutrinosphere to 10 9 cm in order to ease the time 
step constraint. The gravity potential of the central object is 
thereby corrected for the gas mass in the larger excised re- 
gion, and free-outflow conditions are applied at the new in- 
ner boundary. When the SN blast is followed beyond shock 
breakout for several days, the star is embedded in a wind-like 
medium with a density p cc r~ 2 and a constant temperature 
of T = HOOK, and the grid is shifted outward once again to 
cover radii from R^ = 10 10 cm to R a b up to 10 15 cm. 

During the long-time simulations of the SN dynamics we 
follow approximately the explosive nucleosynthesis by solv- 
ing a small alpha-reactions network, similar to the network 
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described in Kifonidis et al. (2003). It consists of the 13 a- 
nuclei: 4 He, 12 C, 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 40 Ca, 44 Ti, 
48 Cr, 52 Fe, 56 Ni, and an additional tracer nucleus, which is 
produced via the reaction rate for 52 Fe(ar,y) 56 Ni within grid 
cells whose electron fraction Y e is below 0.49. This allows us 
to keep track of element formation in regions with neutron ex- 
cess. The network is solved in grid cells whose temperature is 
between 10 s K and 7 x 10 9 K. We assume that at T > 7 x 10 9 K 
all nuclei are photo-disintegrated to a-particles. Such a com- 
position is consistent with the NSE yields that are produced 
by our network solver in the high-temperature limit, because 
the burning network contains a-particles as the only represen- 
tatives of light nuclear species. A feedback from the network 
composition to the EoS and thus to the hydrodynamic evolu- 
tion is neglected. 

2.1. Proto-neutron star core model 

The cooling of the PNS core is described by an analytic 
model that couples the excised core region to the surround- 
ing accretion layer, whose evolution is followed on the com- 
putational grid. We assume the dense PNS core with mass 
M c and radius R c to be approximately homogeneous and its 
EoS to be an ideal T law, P = (T - V)e (with P being 
the pressure, e the internal energy density). Combining the 
total core energy, E c = E g + E- u and the virial theorem, 
£g + 3(F - l)£j + S — 0, we can replace the integrated in- 
ternal energy, E\, and express E c in terms of the Newtonian 
gravitational energy, E g = -|GM 2 //? C , and the surface term, 
S = -4^/? 3 P s , for pressure P s at R c : 

- 3F ~ 4 F S 



3(T-1) 



3(T-1) 



(1) 



The energetic evolution of the core is given by its loss of neu- 
trino energy and compression work done on its surface as 



dE c ■ 

E c = —^= -L v>c - 4kP s R 2 c R c 
dt 



(2) 



where E c can be computed as time derivative of Eq. ([TJ, L vc is 
the total neutrino luminosity, and the second term results from 
the time derivative of the core volume. Instead of setting P s 
equal to the boundary pressure on the hydro grid, we prefer 
to link it to overall properties of the accretion layer. This pre- 
scription is intended to capture the nature of the core-mantle 
coupling but not to constrain the freedom to tune the param- 
eters of the simple PNS core model. We therefore consider 
hydrostatic equilibrium in terms of the mass coordinate m(r), 
dP/dm = -GM/(4nr 4 ), and linearize both sides to obtain 



, GM c m a 



(3) 



where m acc is the mass of the accretion layer that surrounds 
the PNS core, (>0a numerical factor of order unity, and we 
assumed Pq <s P s for the pressure Pq outside of the accretion 
layer. Moreover, in performing the time derivative of Eq. (Q3 
we assumed M c and F to be constant. Combining Eqs. ([TJ — <[3j 
we arrive at 



3r-4 R c £ 6E. dcc 



(4) 



with S = -t,GM c m^lR z and SE acc /5t = GM c ra«/fi c . 
While the first term on the rhs of Eq. (2) describes the lumi- 
nosity increase due to the deepening of the gravity potential 




20 25 30 
ZAMS Mass [M Q ] 

Fig. 1. — Timescale of 90% of the neutrino-energy loss of the forming com- 
pact remnant as function of the progenitor ZAMS mass. Red histogram bars 
indicate successful explosions, grey ones correspond to cases where BHs 
form without a SN explosion, and the green bar marks the 19.8 M Q progeni- 
tor used for the calibration with SN 1987 A observations (see text). The BH 
formation cases correspond to "cooling times" in excess of 6 s, because the 
compact object in our simulations remains radiating neutrinos even when its 
mass nominally exceeds the BH formation limit. This implies that our mod- 
eling does not invoke any assumption about the equation-of-state dependent 
mass limit for BH formation. 



and surface work in the case of PNS compression, the sec- 
ond term accounts for the higher core pressure (and internal 
energy) needed when the accretion layer grows in mass. 

Equation (@]i is used to prescribe the boundary luminosities, 
L V . ;C = ( Vl L ViC , of neutrinos of all kinds with 6 Ye = 0.20, i Ye = 
0.15, i Vx = 0.1625. (This choice corresp onds to a certain los s 
of lepton number from the PNS core, see S check et al.l (120061) . 
and ensures a reasonable evolution of Y e in the PNS mantle 
and surface layers, but does not have much relevance for the 
dynamical evolution.) 

In our simulations m acc is taken to be the mass between the 
inner grid boundary and a density of po = 10 10 gcnT 3 at ra- 
dius ro, where we define m acc = -4-TTr^vopo (accretion means 
a velocity vo < and m acc > 0). The core radius is assumed to 
contract according to R c (t) = R Cii + (R ci - R c f)/(l + t) n with 
^c,i = and R c f being the initial and final radius, respec- 
tively, and t is measured in seconds. With r = 3, n — 3, and 
£ = 0.6 a choice of R c $ = 6 km allows us to reproduce £g7A 
and Mn; 87a of SN 1987 A for progenitors in the 20 M Q range. 
For the simulations discussed below the 19.8 M progenitor 
serves for the calibration of the PNS core model, but the over- 
all results are similar when neighboring stars or a SN 1987A 
blue supergiant progenitor (IWooslev et al.|[l988l) are used for 
the calibration. 

The chosen parameter values lead to typical PNS neutrino- 
cooling times (tgo for 90% of the total neutrino-energy re- 
lease) of 3.5-5.5 s (Fig. Q]). This is shorter than the ~10s 
of emission inferred from the SN 1987 A neutrino events of 
Kamiokande II. However, this detector reported a 7 s gap af- 
ter 8 events in the first two seconds, an d the last 3 event s 
were very close to the detection threshold dHirata et al.lll987l) . 
It is interesting to note that the neutrino signal in all three 
experiments (Kamiokande II, Irvine-Michigan-Brookhaven, 
and Baksan) is compatible with a PNS cooli ng period (ex- 
ponential cooling timescale) of only 4-5 s dRaffeltl 119961 : 
lLoredo & Lambll2002t IPagliaroli et al.ll2009l) . We emphasize 
that our ID simulations, which ignore postshock convection 
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and yield standard values for the mean energies of radiated 
v e and v e ((e) « 10-17 MeV), naturally require overestimated 
postbounce neutrino fluxes to trigger explosions. Since the 
neutrino emission of our PNS model is roughly compatible 
with the SN 1987 A data, we expect that our calibration by 
SN 1987 A explosion properties reflects trends of neutrino- 
driven explosions whose validity for the investigated progen- 
itors in relative comparison holds beyond the considered ID 
setup. We also stress that our description of the core behav- 
ior accounts for the presence of an accretion layer surround- 
ing the (excised) high-density core. The growth of this layer 
depends on the progenitor-specific stellar structure and cor- 
responding mass-infall rate. Its accretion luminosity, which 
adds to the analytically modeled core-neutrino luminosity, is 
also included in our simulations by means of a simple, grey 
transport approximation. 

The progenitor structure therefore influences the post- 
bounce evolution in different ways. The density profile of the 
outer Fe-core, Si-, and O-layers does not only determine the 
mass accretion rate and thus the accretion luminosity; it also 
governs the shock stagnation radius because high/low accre- 
tion rates damp/favor shock expansion. Moreover, according 
to Eq. our PNS-core model connects the neutrino emis- 
sion of the PNS core to the growth of the accretion mantle of 
the PNS. Core and accretion luminosity conspire in reviving 
and powering the shock wave, for which reason the neutrino 
emission properties of the PNS (core and accretion contribu- 
tions) are crucial for the explosion. It is therefore clear that 
the calibration of the PNS-core parameters by the explosion 
properties of a chosen progenitor model cannot be indepen- 
dent of the structure of this progenitor star. The fact that basic 
aspects of our results exhibit some robustness against varia- 
tions of the calibration model is assuring. However, this does 
not imply that the predicted explosion properties could not 
change if the core structure of the SN 1987 A progenitor was 
considerably different from the ~20 M Q solar metallicity stars 
and the blue supergiant model used in our calibration tests. 
Metallicity effects, interaction in a binary scenario, or rota- 
tion and convective boundary effects could have affected the 
mass-loss history and thus could have altered the growth of 
the stellar core compared to the progenitors models consid- 
ered by us for calibration. 

2.2. Progenitor stars 

Density profiles versus radius and mass for a subset of the 
investigated 101 sola r-metallicity progenitors with iron cores 
dWooslev et al .1120021) are displayed in Fig. [2] The models are 
given in 0.2 M steps between 10.8 M and 28 M and fur- 
ther up to 40 M Q in 1.0 M Q steps. We also include two stars 
with 10.0 M Q and 10.2 M . While the density profiles are la- 
beled with the zero-age-main-sequence (ZAMS) masses, the 
corresponding pre-collapse masses are visible in Fig. [3] The 
images show that mass loss is moderate up to 19.6 M Q (ex- 
cept for the 10.0 M Q case) but grows considerably beyond that 
value. 

The nonmonotonic star-to-star scatter is a consequence of 
the turn-on and -off of convective shells during the late evo- 
lution stages. The scatter can be seen in Fig. |4]for different 
structural parameters as functions of the ZAMS mass: The 
iron-core mass, Mp e , is identified with the deleptonized core, 
i.e., the volume where the electron faction Y e < 0.497; the 
binding energy of the overlying shells, Eb, is evaluated at 
the onset of core collapse, and the enclosed mass at the in- 
ner edge of the oxygen-burning shell, Mq < , is defined by the 
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Fig. 2. — Density profiles vs. radius (top) and vs. enclosed mass (bottom) 
at the onset of core collapse for a selection of models from the considered 
set of solar-metallicity progenitors with iron cores. Solid lines correspond to 
ZAMS masses less than 19 M Q , dashed lines to higher values. Note that the 
stellar shell structure and also the high-density core (p > 10 5 gem -3 ) exhibit 
considerable variations with the progenitor mass (see also Fig.U). 
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Fig. 3. — Stellar masses at the onset of core collapse as function of ZAMS 
mass. Mass loss is particularly strong beyond ~20M Q . The red histogram 
bars indicate successful explosions, the grey shaded ones correspond to cases 
where BHs form without a SN explosion, and the green bar marks the 
19.8 M Q progenitor used for the calibration with SN 1987 A observations (see 
text). 

point where the dimensionless entropy per nucleon reaches a 
value of 4. Final ly, the compactness par ameter £2.5 is defined 
as in Eq. (10) of IP Connor & Ottl d2011l) by the ratio of mass 
M = 2.5 M and radius R(M) that encloses this mass: 



h.5 = 



M/Mq 



(5) 



fl(M)/1000km' 
where we consider M as the baryonic mass. For our solar 
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Fig. 4. — Quantities characterizing the core structure of the progenitors at the onset of collapse vs. ZAMS mass: deleptonized ("iron") core mass (top left), 
binding energy of the matter outside of the iron core in units of 1 B = 1 bethe = 10 51 erg (top right) at the onset of collapse, enclosed mass at the bottom of the 
oxygen-burning shell (bottom left), and compactness parameter £?.5 of the innermost 2.5 M G as defined in Eq. |5J (bottom right). Red, grey, and green histogram 
bars have the same meaning as in Fig. [5] 



metallicity progenitors it makes no noticeable difference 
whether we compute £2.5 at the onset of collapse or at the 
time of bounce (as suggested by |Q' Connor & Ottll201 1|) . in 
contrast to the situation for low-metallicity stars, which expe- 
rience less mass loss during their evolution and which develop 
more compact cores. 

The grey histogram bars mark the cases where we found 
BH formation without a SN explosion. No clear correlation 
can be observed with extrema of any of the parameters that 
characterize the density profile of the progenitor core. (We 
note that the binding energy at bounce instead of the value 
at the onset of collapse, though different in detail, does not 
reveal any better correlation.) While many of the BH forma- 
tion cases are located in a local peak of Mp e , E\» Mo < , and 
£2.5 between ~22 M and ~26 M , there are exceptions in this 
interval. Moreover, there are BH formation cases below and 
above this ZAMS mass window with clusterings around 15- 
16M Q , ~18M , and 21-22M . Also the C+O-core and O- 
shell masses, which exhibit a nearly monotonic increase up 
to a progenitor ZAMS mass of 28 M and scattering around 
a high value beyond, do not show any correlations with BH 
formation cases. 

The clearest indication of possible BH formation seem to be 
high values of £2.5, at least local maxima (as visible between 
15 and 16 M and around 18 M ). But there are exceptions 
like progenitors near 21 M , which have relatively high val- 
ues of ^2.5 and £\,- Despite this fact they explode, correlated 
with possessing lower values of Mq k than BH formation cases 
in the close neighborhood. We therefore conclude that the dy- 
namics leading to a successful SN blast does not depend on a 



single property of the progenitor star but on different ones or 
even a combination of them. 

Although of all quantities displayed in Fig.[4]the compact- 
ness parameter £2.5 in the form of local maxima exhibits the 
tightest correlation with BH formati on cases, the corr e spond - 
ing threshold value considered by |Q' Connor & Ottl (1201 ll) . 
i.e., £2.5 > 0.45, significantly underestimates the BH forma- 
tion probability compared to our results. The progenitor vari- 
ability seen in Fig. |H however, defies the specification of a 
better and reliable number for the BH formation limit. In- 
stead, there is a wide range of £2.5 values between ~0.15 on 
the one side and nearly 0.35 on the other side, for which NS 
and BH formation cases occur. Of course, any such conclu- 
sions can be affected by the ID nature of our simulations. It 
cannot be excluded that multi-dimensional effects will shift 
the £2.5 limit for BH formation to larger values. 

3. RESULTS FOR EXPLOSIONS AND REMNANTS 

The stellar core collapse until ~1 0-1 5 ms after bounce 
is computed with the Prometheus -Vertex code includ- 
ing sophisticated n eutrino transport dRampp & Jankal [20021 
iBuras et ai1l2006al) . At that time the bounce shock has con- 
verted to a slowly expanding accretion shock. We then con- 
tinue the simulations of the subsequent accretion phase and 
the possible explosion with the Prometheus version described 
in Sect. [2] The mapping, excision of the NS core, and ap- 
proximate neutrino treatment does not cause any worrysome 
transients. 

3.1. Explosion properties 
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Fig. 5. — Explosion and remnant properties resulting from our parametrized ID neutrino-driven SN simulations: explosion energy (top left), time of the onset of 
the explosion (top right), baryonic mass of the compact remnant (middle left), total release of gravitational binding energy by the compact remnant in neutrinos 
(middle right), and ejected 56 Ni mass (bottom left) as functions of stellar birth (ZAMS) mass. The lower right panel shows the compact remnant mass vs. the 
enclosed mass at the base of the oxygen-burning shell of the progenitor, where the stars possess an entropy jump of varying size. The green histogram bar 
indicates the 19.8 M Q calibration model (see text). While vertical ticks in some panels mark masses where computed models did not explode, grey histogram 
bars reaching to the upper panel edge and arrows in the bottom right panel signal the formation of a BH containing the whole mass of the progenitor at collapse. 
The only exception here is the 37 M Q star, where the explosion expulses ~3.2M G while 4.5 M e of fallback give birth to a BH with 6.5M e . Blue histogram 
segments indicate fallback masses and orange segments the uncertainties of the 56 Ni ejecta masses. The latter uncertainties are associated with inaccuracies in 
the Y e determination of the neutrino-heated ejecta because of our approximative treatment of neutrino transport. 



Explosions can develop in the case of a favorable inter- 
play of mass accretion rate and neut r ino luminosities (e.g., 
iBurrows & Goshv1[l993t Hankal l200lt IPernandezI 120121) . In 
all successful cases compared to failed explosions of neigh- 
boring progenitors, the mass accretion rate is either lower 
during a long postbounce period or decreases rapidly when 
a composition-shell interface arrives at the shock. Shock 
revival occurs when the neutrino luminosity is still suffi- 
ciently high (and thus neutrino heating strong enough) at this 
time. In a large number of successful and unsuccessful mod- 
els the decreasing mass-accretion rate triggers shock oscil- 
lations, which indic a te the proximity to runaway conditions 
dBuras et al.l I2006bt iMurphv & Burrowsl 120081 iFernandezI 
120121) and whose amplification also leads to large-amplitude 
pulses of the accretion compone nt of the driving neutrino lu- 
minosity (see lBuras et al.ll2006bl) . In some stars the explosion 
is fostered by the Si/O interface reaching the shock relatively 
soon after bounce, either due to its location at a smaller mass 



coordinate or because of higher mass accretion rates at earlier 
times, corresponding to a more compact Si-layer. In this case 
the high accretion luminosity seems to be supportive. (More 
information on the time evolution, dynamics, and the neutrino 
emission of our models will be provided in a separate paper). 
In summary, the destiny of a collapsing star does not hinge on 
a single parameter but depends on the overall structure of the 
stellar core. 

Figure gives an overview of the results of our whole 
model set. All displayed quantities exhibit considerable scat- 
ter even in narrow mass windows, which is a consequence 
of the nonmonotonicities of the progenitor structure. Failed 
explosions with BH formation seem to be possible for pro- 
genitors below 20 M , and successful SNe with NS formation 
are found also between 20 and 40 M©. While below 15 M 
all core collapses produce NSs, the investigated progenitor 
set yields several "islands" with preferred BH creation above 
15 M Q . A discussion how BH formation cases correlate or 
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Fig. 6. — Normalized distribution of (baryonic) remnant masses originating 
from the neutrino-driven explosion models of our investigated set of solar- 
metallicity progenitors, weighted by Salpeter's (1955) stellar birth rate func- 
tion. While the main panel shows the fraction of remnants per 0.25 M Q bin, 
the insert provides an enlarged and refined display of the NS mass range with 
histogram bars giving the remnant fractions per 0.10 M B bin. 

do not correlate with the density structure and charac teris tic 
quantities of the progenitor cores can be found in Sect. 12.21 

The energies of the neutrino-driven explosions do not ex- 
ceed 2 x 10 51 erg and 56 Ni production up to 0.1-0. 15 M can 
be expected. Note that our determination of nickel yields 
is uncertain because the neutrino transport approximation 
does not provide very accurate information of the Y e in the 
neutrino-heated ejecta and in the neutrino-driven wind, where 
a sizable contribution to the ejected nickel may come from 
(orange histogram bars in the bottom left panel of Fig. [5]de- 
note the corresponding production masses of the neutron-rich 
tracer nucleus in the network). The onset of the SN blast 
(roughly defined by the moment the shock passes 500 km) 
varies in time between ~0.2-1.2s after bounce, thus includ- 
ing "early" as well as "late" explosions. Later explosions 
tend to be less energetic, because less mass is available in the 
more dilute gain layer for absorbing energy from neutrinos. 
This simple relation for each progenitor is partly masked by 
the large differences of the mass infall rates in the different 
stars: A dense stellar core allows for a high mass accretion 
rate by the shock and a relatively massive gain layer even at 
later postbounce times. Moreover, the neutrino-driven wind, 
which consists of mass shed off the PNS surface by neutrino 
heating, adds a contribution to the explosion power that does 
not depend on the density structure of the silicon and oxygen 
layers. 

The blast-wave energy of the SN emerges from a combi- 
nation of different positive contributions: (a) The total (in- 
ternal plus kinetic plus gravitational) energy of the neutrino- 
heated postshock matter when it begins expansion away from 
the gain radius after shock revival; (b) the energy released by 
the recombination of nucleons and a-particles to heavy nu- 
clei in the expanding postshock material; (c) energy from neu- 
trino absorption in matter that is accreted through the shock 
and channeled toward the gain radius to be partly ejected 
again during a phase of simultaneous accretion and outflow 
after the onset of the explosion; (d) energy carried outward 
by the neutrino-driven wind from the PNS surface after the 
accretion has ended; (e) energy released from nuclear burn- 
ing in the shock-heated ejecta. Two negative effects reduce 
the energy that finally escapes with the SN ejecta: (f) The 



binding energy of the stellar layers swept up by the accel- 
erating explosion shock; (g) energy extracted from ejection 
by mass fallback to the compact remnant. In our explosion 
models we fully account for points (a), (b), (d), (f), and (g), 
while point (c) was estimated t o be potentially releva nt in 
the multi-dimensional situation (Ma rek & Jankal 120091) . and 
effect (e) typically makes only a small contribution to the 
SN explosion energy (the burning of 0. 1 M of Si and O re- 
leases roughly 10 50 erg). Also term (a) is usually unimpor- 
tant, because simulations show that the matter in the post- 
shock layer is marginally unbound (i.e., its total energy is near 
zero) when shock revival sets in and the gas be gins outward 
expansion behind the accelerating blast wave dScheck et al.l 
l2006t IFerna ndez 2012). The dominant positive contributions 
to the explosion energy in our ID models are therefore given 
by terms (b) and (d). The power of the neutrino-driven wind, 
point (d), can account for ~30-70% of the SN explosion en- 
ergy (with higher relative importance in more energetic ex- 
plosions; see the d etailed analysis and especially Fig. C.5 in 
iScheck et al.|[2006tfl 

3.2. Remnant properties 

The NS baryonic masses are in the range of ~1.2-2Af s . 
For the observable gravitational masses these numbers need 
to be reduced by the mass defect associated with the total 
neutrino-energy loss during NS formation, corresponding to 
AE tot /c 2 . Our model estimates of A£ tot are given in the mid- 
dle right panel of Fig. [5] implying corrections roughly be- 
tween 0.1 1 M Q and 0.23 M for the least massive up to the 
most massive NSs. However, these numbers are approxi- 
mative — though in a reasonable range — , and exact binding 
energies depend on the NS EoS and have to be determined 
from detailed general relativistic NS structure models. The 
smallest BH contains 6.5 M and forms in the 37 M Q star by 
late fallback of 4.5 M Q . Since the progenitor has a mass of 
~9.7 M at collapse, the SN explosion ejects roughly 3.2 M . 
All other BHs originate from failed explosions and absorb the 
total mass of the progenitor at collapse (>8.5M ). Although 
the remnant mass follows a trend of growth with the enclosed 
mass at the bottom of the oxygen-burning shell (Fig. [5}, this 
location is no reliable indicator for the fate of the star because 
some models with relatively small silicon core do not explode. 
However, BH formation clusters in regions where local max- 
ima of the NS mass occur. 

Note that we do not terminate the neutrino emission from 
the excised core when the mass of the compact remnant ex- 
ceeds the BH formation limit. Instead, the accreting cen- 
tral object is considered to remain gravitationally stable, for 
which reason its total energy loss can reach AE tot > 5 x 
10 53 erg (middle right panel of Fig. |5) and its cooling time 
can become tgo > 6 s (Fig. [T). However, none of our suc- 
cessful explosions was triggered by this unphysical, long-time 
release of neutrino energy, because all NSs have a baryonic 
mass of <2 M . This also holds for the compact object that 
is temporarily present (i.e., until massive fallback occurs long 
after the launch of the explosion) in the 37 M progenitor. The 
value of the baryonic mass of <2 M is safely below the grav- 
itational mass limit of 1 .97 ± 0.04 M set by the recent pre- 

4 Note that Scheck et al. (2006) provide an analysis of the energy contribu- 
tions only for the first second of the explosion. Their results, however, have 
more general validity because the positive energy injected to the explosion 
after the first second is in most cases roughly compensated by the negative 
energy of point (f). 
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cision measurement of iDemorest et al.l (120101) . according to 
which objects up to at least this mass do not collapse to BHs. 

A gap that seems to be present in the observed remnant 
distribution between NS and BH masses (e.g., lOzeT et al. 
120 10t IValentim et all 120111: iFarr et al.l 1201 lb is compatible 
with our results. This is visible in Fig. [6] which displays 
the normalized (baryonic) remnant mass distribution originat- 
ing from o ur neu trino-driven explosion models, weighted by 
Salpeter's (119551) stellar birth rate function (with power-law 
index -2.35). We obtain a relatively wide distribution of NS 
masses (see insert). Note that these are birth masses result- 
ing from single, nonrotating (or at most slowly rotating) stars, 
for which effects like mass transfer and mass loss associated 
with possible binary interaction during the progenitor evolu- 
tion have been ignored. Accretion of the NSs in binaries will, 
of course, also modify this NS birth-mass distribution. 

The fallback mass that adds to the initial NS mass is shown 
in blue in the middle left panel of Fig. [5] Fallback is bigger 
(up to ~0.2M Q ) for the lower-mass progenitors where an ex- 
tended hydrogen envelope leads to a significant deceleration 
of the explosion shock and a stronger reverse shock. The dis- 
tribution of Fig.|6]implies that 22.5% of the solar-metallicity 
progenitors produce BHs at the end of their lives. We point 
out that in the considered sample of 101 solar-metallicity pro- 
genitors only the 37 M star produces a fallback SN with de- 
layed BH formation. In this case the black hole swallows a 
baryon mass of ~6.5M Q while the explosion ejects the re- 
maining ~3.2M of the pre-collapse mass of the star. 

Compared to the neighboring, exploding stars of 36 M Q 
and 38 M , the 37 M progenitor possesses a higher enclosed 
mass m(r) and higher gravitational binding energy E\,(r) in 
the central 15,000 km. This is reflected by its larger values 
of E b (M Fe ), M 0< , and £ 2 .5 in Fig. |4j although M Fe is slightly 
smaller. For this reason the 37 M Q model explodes only very 
late (f exp ~ 1 s) and with low explosion energy (Fig. |5). Be- 
cause of the late onset of the blast-wave acceleration the neu- 
trino energy injected to the explosion and the shock velocity 
become considerably lower than in the 36 and 38 M Q cases, 
in fact too low to unbind the major part of the star. It is inter- 
esting to note that the 39 M progenitor is even more extreme 
in all of its characteristic parameters displayed in Fig. |4] Cor- 
respondingly, it does not achieve to explode and collapses to 
a BH as a whole. In contrast, the 34 M model is also locally 
extreme in some of its structural parameters, but its compact- 
ness £2.5 is lower and it explodes much earlier than the 37 M 
star. The latter therefore defines an intermediate, seemingly 
rare situation in which the combination of a high binding en- 
ergy and core compactness of the star and a late onset of the 
explosion and low explosion energy favors the fallback of a 
large fraction of the progenitor mass. 

4. CONCLUSIONS 

We performed simulations of neutrino-powered explosions 
in spherical symmetry for a large set of solar-metallicity, 
nonrotating stars with Fe-cores. We calibrated an analytic 
neutrino-cooling model of the PNS core such that the ob- 
served explosion properties of SN 1987 A were reproduced 
for progenitors around 20 M Q . This yields numbers for the to- 
tal energy release and cooling timescale of the newly formed 
NS in rough agreement with the SN 1987 A neutrino mea- 
surement, thus confirming the overall consistency of our ap- 
proach. 

The neutrino-driven explosions turned out to be fostered by 
"steps" in the stellar density and entropy profiles, which re- 



duce the mass-infall rate and ram pressure of the accretion 
flow and thus allow the shock to expand. Explosions set in 
either by a continuous acceleration of the shock or after a 
phase of shock oscillations with growing amplitude. The on- 
set times range from only 0.1s to more than 1 s after core 
bounce, meaning that depending on the progenitor the explo- 
sion can occur relatively early or with a considerable delay. 
This conflicts with the claim that neutrino-driven explosions 
should generally be launched soon after bounce in order to 
account for the obser ved gap between highe st NS masses and 
smallest BH masses dBelczvnski et alj|201 1[) . 

Reflecting the "erratic" structure differences of the progen- 
itor models, the explosion properties and remnant masses ex- 
hibit highly nonmonotonic variability even in narrow ranges 
of the ZAMS mass. This might explain the considerable 
differences of the explosion energies and nickel masses in- 
ferred for observed SNe from lightcurve and spectral analy- 
ses as function of the ZAMS mass of the progenitors (e.g. , 
lUtrobin & Chugall l20TTt iSmartt et al.l 120091: iSmarttl 120091) 
without invoking additional pre-collapse degrees of freedom 
like rapid core rotation. However, the large variations of ex- 
plosion and remnant properties even for small changes of the 
progenitor ZAMS mass should be taken with caution. The 
mixing length or mixing length plus overshooting descrip- 
tion of convection used in the employed progenitor models 
is extremely unstable to numerical effects in the late burning 
stages. While the trends we found in dependence of iron- 
core mass, core-binding energy, and compactness are likely 
to be valid, the large structural differences of progenitor mod- 
els with similar ZAMS masses probably do not reflect phys- 
ical reality. Instead, late stage convection is subject to sub- 
stantial dynamic instability, for which reason there may be 
stochastic variations in the core properties at any given mass. 
An exploration of this possibility, however, has to be deferred 
until multi-dimensional progenitor models, self-consistently 
evolved to the onset of core collapse, become available. 

Neutrino heating seems hardly able to provide explosion 
energies significantly higher than 2 x 10 51 erg and 56 Ni yields 
in excess of 0. 1-0.2 M Q . Observed SNe with higher explo- 
sion energies and nickel production (e.g.. lTanaka et al.ll 19991) 
are therefore likely to be driven by a mechanism that is differ- 
ent from neutrino heating and could involve magnetohydro- 
dynamic processes. 

The distribution of NS birth masses (including fallback) 
spans from ~1.2M to ~2.OM in baryonic mass with a 
broad maximum between 1.4 M and 1.7 M . BH forma- 
tion occurs in certain mass intervals and seems to be pos- 
sible not only for progenitors with ZAMS masses beyond 
20 M Q but also between 15 M and 20 M Q . With a single 
exception (the 37 M progenitor) all BHs swallow the pro- 
genitor completely. This suggests that fallback SNe could be 
rare cases in solar-metallicity environments although roughly 
23% of the progenitor population end as BHs. Because of 
the large mass loss of stars with Mzams > 15 M Q , the fall- 
back is strongest (up to ~O.2M ) for stars below ~2OM . 
Correspondingly, we find a wide gap between the maxi- 
mum NS birth mass (around 1.8 M gravitational mass) and 
the minimum BH mass (roughly 6M Q ), compatible with 
conclusions from observed compact remnant masses (e.g., 
lOzel et al.l l20Tot iValentimetaLl |20Tlt IFarr et al.l l20TTb . It 
should, however, be noted that our predictions do not account 
for binary effects prior and after the SN explosion, nor do we 
include NSs formed from ONeMg-core progenitors. A direct 
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comparison with observations therefore requires caution. 

The compactness parameter £2.5 introduced by 
lO' Connor & Otj d201 ll) provides the best indicative quantity 
based on the pre-collapse structure for the final fate of a star, 
i.e., whether it is likely to end as NS or BH. However, our 
models fail to explode significantly below the discriminatin g 
value of £2.5 = 0.45 considered by |Q' Connor &~Ottl (120111) . 
If this result were valid more generally, in particular also 
for lower-metallicity progenitors, it will h ave implications 
for the discussion of G RB progenitors dO' Connor & Ottl 
1201 It iDessart et aLll2012l) . Instead of depending on an exact 
bifurcation value of £2.5 we see BH formation occurring pref- 
erentially (but not exclusively) in regions with local maxima 
of £2.5 and NS as well as BH formation to be possible for 
£2.5 between ~0.15 and ~0.35. The success or failure of the 
explosion mechanism obviously depends on the progenitor 
structure (determining the time-dependent mass-accretion 
rate, shock radius, and accretion luminosity) in a complex 
way and cannot be predicted exactly on grounds of a single 
parameter. 

So far we have investigated only solar metallicity pro- 
genitors from one modeling group, but we plan to ex- 
tend our studies also on sets of progenitors from other 
groups and for different metallicities. Naturally, the spher- 
ical symmetry of our simulations is a major constraint, 
since the explosion mechanism has b een recognized to 
be intrinsically mu lti-dimensional (e.g., iJanka et ail 120071: 
iBurrows et al.ll2Q07h . Nonradial fluid flows in the gain layer 
have a supportive influence and thus r educe the critical lumi- 
nosity for driving the explosion (e.g.. IJanka & Mullerlll996t 
Murphy & Burrowsll2008t iNordhaus et al.l l2010t lH anke et al.l 
201 U iMurphv & Meakinl hoi ll) . With a higher neutrino- 



heating efficiency lower values of the early postbounce lumi- 
nosities from our analytic PNS core-cooling model would be 
needed to reproduce the explosion properties of SN 1987 A. 
This, in turn, would stretch the PNS cooling timescale, which 
in our simulations tends to be on the lower side of the range 
compatible with the SN 1987 A neutrino data. 

Nevertheless, this would only mean a recalibration of the 
time decay of the core-neutrino source adopted for our mod- 
eling, but it is unlikely to fundamentally alter qualitative as- 
pects of the behavior of the different progenitors in relative 
comparison to each other. Therefore we are hopeful that 
our basic conclusions will retain their validity also in the 
multi-dimensional case even if quantitative changes (for ex- 
ample lower explosion energies by less powerful neutrino- 
driven winds in low-mass progenitors) could be the conse- 
quence. It also cannot be excluded that some of the extreme 
sensitivity of ID explosions to details of the progenitor struc- 
ture might be moderated or disappear in multi-dimensional 
models. It is interesting to note, however, that the relatively 
early explosions of the 11.2M Q and 27 M progenitors of 
our model set and the later explosion of progenitors around 
15M are fully compatible with recent res ults of 2D mod- 
els with sophisticated neu trino transport (cf. iMarek & Jankal 
120091: lMuller¥aniMT2allbl) . 

Our findings challenge canonical thinking, e.g . of the pro- 
genitor masses leading to BH formation (e.g., IZhang etafl 
120081: lO'Connor&Otj 1201 ll) . and they raise doubts about 
claims in the literature that the possible gap between NS and 
BH masses re quires explosions that set in relatively early after 
core bounce dBelczvnski et all 1201 U iFrver et al.ll2012l) . We 
have shown that a highly variable onset time of the explosion 



and little fallback for massive progenitors (i.e., a low occur- 
rence rate of fallback SNe) naturally leads to the mass gap. 
On the other hand, a long delay of the beginning of the ex- 
plosion is questioned by the recent Bayesian an alysis of the 
NS m ass distribution in double NS systems by iPejcha et al.l 
(120121) . These authors conclude that for solar metallicity stars 
the explosion most likely develops near the edge of the iron 
core at a specific entropy of S /(A^a^b) ~ 2.8. It is hard to 
understand how the explosion mechanism could dynamically 
couple to this specific pre-collapse value of the entropy, be- 
cause the core density and entropy profiles are highly variable 
for the different progenitors, they are modified by the nuclear 
burning during the collapse phase, and the neutrino-driven 
mechanism in spherical symmetry as well as in the multi- 
dimensional case is sensitive to the mass-infall rate rather 
than the pre-collapse entropy structure. For these reasons 
the mass-accretion rate and the accretion luminosity exhibit 
large differences between different progenitors when com- 
pared for specific values of the entropy, for which reason a 
wide range of explosion times ap pears plausible. Interest- 
ingly, however. IPejcha et al.l(l2012l) favor a NS formation sce- 
nario with no or very little mass fallback, which is consistent 
with the tendency of our models compared to piston-driven 
explosion models . Correspondingly, the probability distribu- 
tion computed bv lPeicha et alj d2012l) with our predicted rem- 
nant masses has little power for NS masses above ~1.6M 
even when our fallback masses are included. Nevertheless, 
the theoretical distribution is shifted by ~0. 1 M to higher 
masses relative to the observational double NS data. If this 
mismatch cannot be attributed to an overestimated delay time 
of the explosion in our supernova models as argued above, it 
might be connected to the low-number statistics of the empir- 
ical data or might point to smaller iron-cores than present in 
the solar-metallicity progenitors of our study. 

Certainly, our simulations can mean only a very first step 
of systematically exploring the consequences of neutrino- 
powered explosions for a wide variety of progenitors. Nev- 
ertheless, our results already shed new light on a number of 
important astrophysical questions connected to how the prop- 
erties of SN explosions and their remnants connect to those 
of the progenitor stars. While the fundamental dependences 
of explosion and remnant properties on the progenitor struc- 
ture are likely to possess more general validity, the detailed 
functional variation with the progenitor mass is subject to a 
number of caveats. These caveats are, in particular, the em- 
ployed set of progenitor models with their large structural 
variations even for similar ZAMS masses, the calibration of 
the neutrino-cooling models of the PNS core based on a cer- 
tain choice of the SN 1987 A progenitor model, and the con- 
straint of the explosion modeling to spherical symmetry. Fu- 
ture work will have to explore the consequences of these con- 
straints and deficiencies. 
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